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ABSTRACT 


A method for recursively computing the coefficient 
of the time-dependent f and g series solution to 
the two-body problem is derived. To find the coeffi- 
cients of the terms in the f and g series, time 

-z 

derivatives are taken of the function h = p./r . I* 1 

this formulation another variable 0 = hr is intro- 
duced, and a recursive relation for the derivatives of 
h is found in terms of the derivatives of the radial 
magnitude and of the preceding derivatives of h. A 
recursive relation for the derivatives of r is then 
derived. After the time derivatives of h are deter- 
mined, the coefficients in the f and g series are 
found from well-known formulas. The convergence prop- 
erties of the f and g series are studied using up 
to 102 terms. The time interval of convergence is com- 
pared with the computed time radius of convergence of 
the series with excellent agreement. 
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RECURSIVE COMPUTATION OF THE COEFFICIENTS OF THE TIME -DEPENDENT 

f AND g SERIES SOLUTION OF KEPLER IAN MOTION AND A 

STUDY OF THE CONVERGENCE PROPERTIES OF THE SOLUTION 

By Victor R. Bond 
Manned Spacecraft Center 


SUMMARY 


This report formulates a method for recursively computing the coefficients 
of the time -dependent f and g series solution to the two-body problem. To 
find the coefficients of the terms in the f and g time series, time deriva- 
tives are taken of the function, h * \x /r^. This ordinarily becomes algebrai- 
cally prohibitive after only a few derivatives. In the formulation presented 
in this paper another variable, 0 = hr, is introduced, which gives a recursive 
relation for the derivatives of h in terms of the derivatives of the radial 
magnitude, r, and of preceding derivatives of h. Recursive relations for the 

derivatives of 0 and r are then derived. After the time derivatives of h 

are determined, the coefficients in the f and g series are found from stand- 
ard formulas. 

The data presented provide a comparison of the non-dimensional angular 
momentum fg - gf computed by this recursive method with the theoretical 

value of unity for the angular momentum. This comparison indicates the number 

of terms which should be used in the f and g series for them to represent 
accurately a solution to the two-body problem for a given set of initial con- 
ditions. 

The results show that there is a point beyond which the use of additional 
terms in the series does not extend the time interval over which the series 
solution is valid. This time interval may be computed from the formula for the 
radius of convergence of the f and g series. The regions of convergence 
that are apparent from the results are in excellent agreement with the computed 
values of the radius of convergence. It is demonstrated by four examples that 
the use of 30 terms in the series maintains a high degree of accuracy for a 
time interval of approximately one-half the computed radius of convergence. 


INTRODUCTION 


The solution to the two-body problem, or Keplerian motion, may be expressed 
explicitly as a function of time through the use of the f and g series. As 


stated in reference 1, this solution has "been known for some time, having first 
been used by Lagrange in 1869. The problem encountered in using this solution 
is that the coefficients in the series have been very difficult to obtain except 
for about the first six terms. In recent literature (refs. 2 and 3) authors 
comment that six terms are sufficient for short time intervals and that the 
higher terms are tedious or impractical to obtain. Reference 4 indicates that 
a general or recursive formula for generating the coefficients exists but is 
obscure; no attempt is made to present the recursive relations. 

Reference 1 describes a digital computer program which generates up to 
27 coefficients for use in the f and g series. This program uses a re- 
cursive formulation first discovered by Cipolletti in 1872. This formulation 
is recursive in only a limited sense because the coefficients of a given order 
depend upon the derivatives of the coefficients of the previous order. The 
coefficients beyond the sixth or seventh are still impractical to obtain by 
hand computation because there is no recursive formulation for the derivatives 
of the coefficients in Cipolletti’ s method. 

This report includes the derivation of a recursive method for generating 
the coefficients which does not require derivative operations. It also dis- 
cusses some of the convergence properties of the f and g series and pre- 
sents several numerical examples. 


SYMBOIB 


a semimajor axis 

a, coefficients in equation for f(t) 

it 

b, coefficients in equation for g(t) 

K 

c^ coefficients in equation for h(t) 

E eccentric anomaly 

e eccentricity 

f(t) series defined by equation (2) 

g(t) series defined by equation (3) 

H the hyperbolic equivalent of E 

h(t) series defined by equation (8) 

i index inferred from equations (l6) and (17); \J -1 in appendix A 

j index inferred from equations (4) and (5) 
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L angular momentum per unit mass 

/ • 

L the non-dimensional angular momentum, fg - gf 

In natural logarithm 

M mean anomaly 

N the hyperbolic equivalent of M 

n index inferred from equations (10), (l6), and ( 17 ) 

n maximum number of derivatives of h desired 

max 

R-d radius of attracting body 

JD 

r position vector 

r magnitude of r 

T time of peri-apsis passage 

t time 


v velocity vector 

a parameter defined by equation (A7) 

M» universal gravitational constant times mass of attracting body 

p radius of convergence 

radius of convergence for hyperbolic orbits 
0 variable defined by equation (12) 

Other Notations : 

( n^* 1 derivative of ( ) with respect to time 

( ) derivative of ( ) with respect to time 

( ) ( ) evaluated at initial time t : or indicates the first coefficient 

x o o 

in the series defined "by equations ( 2 ), (3)> and (8) 




the binomial coefficients 


( )* value of ( ) at first singular point of the f and g series 
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DERIVATION 


This formulation of the f and g time series solution to the two-body 
problem will follow to a large extent that of Brouwer and Clemence (ref. 5), 
except that vector notation will be used for convenience. 


Tae solution to the two -body problem may be written 

r (t ) = r Q f(t) + v q g(t) 

r and v are the initial position and velocity vectors. The functions 
— o — o 

f(t) and g(t) are given by 

CO 

k 


(1 ^ 


f(t) 'Z a X (* ' *°) 

k=o 

CO 

g(t) -Z (* - t o ) 1 

k=:Q 


The a, and b n are obtained from the relations 
k k 


a 


- 1 


j+2 " U + D( J + 2) 


u 

Z 

k=o 


C k a j-k 


( 2 ) 


(3) 


(4) 


b. 


- 1 


j+2 ~ (j + l)(j + 2) 


C k b j-k 


k=o 


for j ^ o. The coefficients a , a , b . b are 


a = 1 
o 

a^ = 0 

b « 0 
o 

h = 1 




j 


(5) 


(6) 


The c are derivable from the Taylor's expansion for the term p/r^ 


which appears in the equations of motion 


b 



(7) 


r + r = 0 


Define the function 


h(t) = ^ 


( 8 ) 


by expanding h(t) about the epoch t 


h(t) C[i (t . t o ) n 

n=o ' / 


(9) 


where 


c = 
n 


h <V 

nl 


(n) 


( 10 ) 


The subscript zero will now be dropped from the notation, since all co- 
efficients and derivatives will be assumed to be evaluated at t . 

o 

Recursive Formula for Time Derivatives of h 

The time derivatives of h must now be obtained in order to develop the 
solution. 

The first derivative of h is 


h (D , ^ T (i) , . r (D 


(ii) 


At this point introduce the variable 


<f = hr 


(i) 


Equation (ll) becomes 


h 


(i) . _ 


Take the second derivative of h from (13) to obtain 

h (2) = _ 5 r^ 1 ) . ft/ 1 ) 


( 12 ) 


(13) 


( 14 ) 
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Now use (13) to eliminate in (l4) 


(15) 


In successive derivatives a similar elimination process is followed to obtain 

h ^ + 2r^ TcS 2 ^ + r^ h^) 


The numerical coefficients of the terms in these expressions may be found from 
the binomial coefficient relation. For the nth derivative 


h 


(n) 


-1 

r 


30 


(n-l) 


n-l 


n-l 

i 


r (i) h (n-i)| 


i=l 


( 16 ) 


Recursive Formula for Time Derivatives of 0 

The derivative 0^ n must now be found in order to evaluate properly 
( 16 ). Applying Leibnitz's rule for the differentiation of a product to equa- 
tion ( 12 ), the nth derivative of 0 becomes 


0^ n ) = hr^ n+1 ^ h^ r (n-i+l) 

i=l 


(17) 


Equation ( 17 ) may now be used to eliminate ^ 

stitute n-l for n in ( 17 ) and obtain 


in equation (l 6). Sub- 


0 (n - l} = hr (n) +£ ( n “ 1 ) h (i) r (n - i} 
i=l 


( 18 ) 


Now substitute ( 18 ) into (l 6 ) and obtain 


h 


(n) 



n-l 




(n-i) 

r 



(19) 
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Recursive Formula for Time Derivatives of r 

The derivatives of the radial magnitude are now needed to complete the 
form ula tion. The first derivative of r is found from the initial conditions 


.(1) ~Q * ~Q 


The second derivative is found from the equation of motion 


.(2) w _ Ul. 


3 2 
r^ r 


where 


L = r Xv 
-o — o 


Equation (21 ) may he written by using (8) 


.(2) _ £ h _ LL_ 
r 2 


The third derivative of r is 


r (3> . £ h (D + auli! 

n ~3 


or by using (8) and (12) 


r (3> . £ h (l) + 2 „ 

11 

The higher derivatives of r follow directly 

r^ = — h^ + 20^ 

M. 


The r ( n+ 2 ) derivative of r may now be written 

r (»+2) , I ? h (n) + 20 (n-l) 

u 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


(24) 


7 



(25) 


or by using equation (l8) 

r («+2> , £. h (n) + 2 

M- 

for n ^ 1. ■ 


( n - 1 ) h (l) r ^- 11 + 


1=1 


Method of Computation 

Given: r,v,u,t-t,n 

-nr — o ^ o max 


(a) 

Compute h, r "*4 and from equations (8), 

0>) 

Set n = 

1 



Compute 

hW 

from ( 19 ) 


Compute 

r (3) 

from (25). 

(c) 

Set n = 

2 



Compute 

h< 2 > 

from ( 19 ) 


Compute 

rC) 

from (25 ) 

(a) 

Continue 

this 

scheme until the desired number 


max 


tives of h are found. 

(e) Compute the c , n = 0, to n 


from equat ion ( 10 ) . 


(f ) Compute the coefficients a and b , k = 0 to n + 2 from 

k k max 

equations (4) and (5). 

(g) Compute f(t) and g(t) from equations (2) and (3). 

(h) Compute the position vector r(t) from (l) and also the velocity 
vector from the time derivative of (l). 

This completes the solution, 

CONVERGENCE OF THE SERIES 


In reference 6 Sconzo and Hale shoved that the radius of convergence, that 
is, the time between epoch time t and the time of the nearest singular point 

of f and g, could be found from the expression 
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p = 


— < M 


In 


i + s/i 


- \ a 7 


1 

.n2>2 


(26) 


where 


0 < e < 1. -it ^ M ^ jt, and a > 0 
9 o 9 


This relation is obviously limited to elliptical orbits, since for hyperbolic 
orbits, a < 0. 

The semimajor axis a, eccentricity e, and the mean anomaly M q may 
be computed from the formulas 


1 = 2_ 
a ~ r 


v . v 
— o — o 

M. 


e cos E = 1 - — 
o a 


r . v 
•n — o — o 
e sxn E = — — 

0 \A±& 


M = E - e sin 
o o 


111 E ° = 751 ° - T ) 


( 27 ) 

(28) 

( 29 ) 

( 30 ) 


The implication of equation (26) is that the convergence of the series de- 
pends not only upon the shape of the orbit as defined by e and a, but also 
upon the initial time as measured from the time of periapsis passage T. 


The radius of convergence for hyperbolic orbits is given in appendix A. 


M3MERICAL EXAMPLES 


Four examples are presented to illustrate some of the results of the re- 
cursive formulation. In each example, the values of the quantity L = fg - gf 
are computed several times, each time with a different number of terms ranging 
from 6 terms to 102 terms. The quantity L is a constant of the motion as 
shown in appendix B and should have a constant value of unity. Any deviation 
of L from unity must result from either truncating the series before a suf- 
ficient number of terms have been taken, or from extending the time interval 
beyond the radius of convergence of the series. 
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Example 1, low altitude nearly circular orbit about the moon: 
r^ = (-1012. 4570, -51.263872, -20.120039) n. mi. 

^ = (-287.96060, 4914.2673, 1967.3377) ft/sec. 

Example 2, at insertion into translunar orbit from earth parking orbit: 
r^ = (3091.8028, 1633.2175, 883.53^7) n. mi. 

* (-14646.307, 27051.882, 17921.139) ft/sec. 

Example 3, approximately 46 hours from insertion into a translunar trajectory: 
r = (-161265.14, -20351.149, -5044.6929) n. mi. 

— o 

v * (-3132.8173, -1023.4259, -501.70741) ft/sec. 

— o 

Example 4, at lunar sphere of influence on a translunar trajectory, going 
toward the moon: 

= (25135.706, -20187.383, -11280.829) n. mi. 
v = (-3090.2697, 2125.8987, 1198.1489) ft/sec 


In examples 2 and 3, and v are referenced to the geocentric system 

and in examples 1 and 4, and are referenced to the selenocentric 

system. 

All four examples, seen in tables I to IV, show that the accuracy of the 
series solution improves as the number of terms in the f and g series in- 
creases for a given time from epoch, as would be expected. It is also seen 
from tables I to IV that the series solution diverges rather sharply near the 
computed values of the radius of convergence. The addition of more terms to 
the series will not increase the validity range of the solution near the com- 
puted value of the radius of convergence. 

An interesting point that can be seen from tables II to IV is that 
i/ may actually be closer to unity with only a few terms than it is with a 
larger number of terms; but this occurs only at values of time greater than the 
computed radius of convergence. The numbers for L beyond the radius of con- 
vergence have little if any meaning. 

The radius of convergence for example 1 was found from equation (26) to be 
P =3.69 hours. The eccentricity of this orbit is 0. 00001 O32, so that it is 
very nearly circular. For a purely circular orbit where the eccentricity is 
exactly zero, the radius of convergence would approach infinite time. 
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The radii of convergence p, for examples 2 and 3 ("both are examples of 
nearly parabolic trajectories) , were computed from equation (26). For example 2 
it was found that p = 0. 2280 hour, and for example 3 it was found that 
p - 45.71 hours. These times are consistent with the results shown in tables II 
and III. For example 2 the epoch chosen was 0.03941 hour after perigee, and, 
for example 3 the epoch chosen was 45.71 hours after perigee. The contrast of 
these two examples points out the implication of equation (26), that is, that 
the radius of convergence depends upon the relation of epoch to the time of 
periapsis passage. The radius of convergence increases as t Q - T, or M q 
increases. 

Example 4 is a case representing a hyperbolic orbit. The epoch chosen was 
13.51 hours before pericynthion passage, and the radius of convergence was 
13.51 hours. 

The number of terms required for the solution is of course dictated by the 
accuracy requirements of the problem to be solved. It should be noted from 
tables I to TV that for 30 terms, eight-digit accuracy is maintained for a time 
equal to about one-half of the computed radius of convergence. The time inter- 
val over which the series converges increases with the addition of more terms, 
but the change in the time interval decreases. This means, of course, that the 
addition of more terms will increase accuracy, but not always by a significant 
amount. For example, table I shows that eight-digit accuracy is maintained to 
a maximum of 2.20 hours by using 54 terms or by using up to 102 terms. The 
accuracy that is obtained, of course, depends to a great extent on the accuracy 
of the computer that is used in the computation. The numerical results pre- 
sented here were obtained on a 10-digit electronic computer. 


All four examples were computed by using the non-dimensional ized method 
such that the initial position was in units of the radius FL, of the attract- 

-D 

ing body; the initial velocity in units of the circular satellite velocity 




and time in units of 



1/2 


This allows the gravitational pa- 


rameter [I 


to be set equal to unity in the equations. 


CONCLUDING REMARKS 


A recursive scheme for computing the coefficients in the f and g series 
solution to the two-body problem has been presented. Four examples demonstrate 
the time interval over which the series solution is valid from 6 terms up to 
102 terms. This time interval is compared with the radius of convergence of 
the f and g series, and excellent agreement is obtained. It is demonstrated 
by four examples that a high degree of accuracy is maintained for a time inter- 
val of approximately one-half the computed radius of convergence by using 
30 terms in the series. 
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APPENDIX A 


RADIUS OF CONVERGENCE FOR HYPERBOLIC CASES 


In reference it was shown that the radius of convergence of the f and 
g series for the elliptical case could he expressed by 



(Al) 


where M* is the value of the mean anomaly at the first singular point of the 
f and g series. By analogy, the radius of convergence for the hyperbolic 
case may he derived by expressing equation ( Al ) as 



(A2) 


where N is the hyperbolic equivalent of the elliptical mean anomaly M, and 
is given by 

N = e sinh H - H (A3) 


The value of N at the first singular point of the f and g series is N*, 
and the value of N at epoch time is JSf . The hyperbolic equivalent of the 

elliptical eccentric anomaly E is H. 

The singularities of equation (A 3 ) are given by 

^ = e cosh H - 1 = 0 (A4) 


This follows since the singularities of the function H(N) are given by 


<15 - 

dN 


(A5) 


Since for the hyperbolic case 
solutions of equation (Ak) for H* 


e > 1 , and also since cosh H ^ 1 , the 
can only be in imaginary. Setting 


H* = ia 


(a6 ) 


•where 


i 


- /u 
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I 


and substituting this expression into (a4), 

e cosh (ia) -1=0 

or, since cosh (ia) = cos a 


cos a = — 
e 


(A7) 


Substituting equation (A 6) and equation (A7) into equation (A3), the value 
of F* at the first singularity is 

F* = — — — sinh (ia) - ia 
cos a v 

But since sinh (ia) = i sin a 


F* = i(tan a - a) 

Substituting equation (A8) into equation (A2) yields 


(A8) 


P H / n 


i(tan a - a) - N 


The absolute value of the complex variable p is given by 

n 


P H = 


-a" 


N + (tan a - a) 
o 


^l 2 


(A9) 


The semimajor axis a may be computed as in the elliptical case from 
equation (27)- The eccentricity and are computed from 

r 

e cosh H = 1 - — (A10) 

o a 


r . v 

e sinh H = .° (All) 

° V-Ma 


B 0 = e sinh « o - H o - - T) (A12) 

u a 

The parameter a is computed from equation (A7). 
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APPENDIX B 


THE NON-DIMENSIONAL ANGULAR MOMENTUM, fg - gf 


The quantity fg - gf may he shown to he a constant of the motion hy 
taking the cross product of equation (l)'with its time derivative 


v(t) = + 


(Bl) 


Taking the cross product, 

r X v = (r Q f + v Q g) x (r Q f + Z Q g) 

= X ^o) (fg - gt) 

Since the angular momentum r Q x X a mus ’t remain constant, 

r X v = r X v 


Therefore, 


’g - gf = 1 


(B2) 


It is seen also that the ratio 


r X v 


r X v 
— o — o 


= fg - gf 


(B3) 


This quantity is defined hy L and is called the non-dimensional angular 
momentum. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, September 9 , 19^5 
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TABLE I.- NON-DIMENSIONAL ANGULAR MOMENTUM AS A FUNCTION OF TIME (EXAMPLE l) FOR A LOW ALTITUDE NEARLY CIRCULAR ORBIT ABOUT THE MOON 


£a = IOI3.93 n. mi. ; e = 0 . 00001032 ; t Q - T = 0.8017 hr; p = 3-694 hr] 


Time, hr 

S ’L / (6 terms) 

L / (l 8 terms) 

L # (30 terms) 

L # (42 terms) 

L ; (54 terms) 

L ; (66 terms) 

L f (78 terms) 

L / (90 terms) 

L # (102 terms) 

0,00 

1 .'inno'iooo 

i.uo iooooo 

1 .onoooono 

1 .00000000 

1.00000000 

1 .00(10 0000 

1 .noOOOOOO 

1 .nooooooo 

1 .00000000 

o.ao 


1 .oooooono 

1 .onoooono 

1 .oooonoon 

l .00000000 

1 .00000000 

1 .nooooooo 

1 .nooooooo 

1.00000000 

0,40 

n.9H1 82 T-i 9 

.i.uo>'onono 

1 .nooooono 

1 .00000000 

1 .00000000 

1. nooooooo 

1 .nooooooo 

1 ,00000000 

T. nooooooo 

0.60 

0.82045322 

l.tMHonono 

1.00000000 

1 .nooooono 

1 .00007000 

1.00000000 

1 .nooooooo 

1 .oonooooo 

1.00000000 

0.80 

n. 207203)2 

I. 0 O 0 OOOOO 

1 .onoooono 

) .nooooono 

1 .nooooooo 

1 , nooooooo 

1 .00000000 

1 ,00000000 

1.00000000 

1.00 

-0,96464*. >8 

<i *99999967 

1.00000000 

^ .oonooooo 

1 .00000000 

1.00000000 

1.00000000 

1 .oonooooo 

i .00000000 

1.20 

-0.99H97217 

0.9999916b 

1 .onoooono 

1.00000000 

1 .00000000 

1 .00000000 

1 .00000000 

1 ,00000000 

i .nooooooo 

1.40 

7.40513006 

U. 99990544 

1 .nooooono 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1,00000000 

1 .nooooooo 

1.60 

4S.OH339 79R 

U .99951540 

0.99999999 

1.00000000 

1.00000000 

1.00000000 

1 .nooooooo 

1,00000000 

i. nooooooo 

1.80 

l bn . 455967 1 3 

1 * u0'?6474 

0.99999865 

1.00000000 

1.00000000 

1.00000000 

1 .00000000 

1 ,00000000 

i .00000000 

2.00 

44n »^71 

l .u? 14741? 

0.99996624 

1 .00000003 

1 .00000000 

1.00000000 

1.00000000 

1,00000000 

1 . nooooooo 

2.20 

1 Ob9. 925b 7 t>*9 

1 . 1 b 484599 

0.9994904) 

1.00000072 

1 .oooonnoo 

l .nooooooo 

l.oooonooo 

1.00000000 

1.00000000 

2.40 

23H3.41 4M JR/9 

1 .7 0*419094 

0.99590675 

1.00000269 

1 .OOOOOOU 

1 .00000001 

1*. 00000001 

1.00000001 

1 .00000001 

2.60 

4632.25599**1 

2. 77^97663 

0.99267234 

0 .999*7866 

1.00001480 

0.99999974 

0.99999997 

0.99999997 

0.99999997 

2.80 

8757.63835345 

1 .05990198 

1.23112709 

«. 98096202 

1.00079465 

0.99997926 

i.oooonois 

0.99999999 

0.99999997 

3.00 

15737.5056*864 

-19.16049293 

4.20531117 

*‘.65112437 

1.02244277 

0.99945287 

0. 99992753 

1 .00001273 

0.99999881 

3.20 

27098.9654*9 46 

-149.92193010 

24.85781020 

-2 .65 7 879 1 2 

1.22616860 

1.05622703 

0.97741093 

1.00469948 

0.99932236 

3.40 

44990.32841 1 ID 

-92 1 .60817027 

1 ?9. 45372920 

-9.72631443 

-S. 43368079 

7.5112 U 066 

-1 .84211140 

1 .84966776 

0.85383668 

3.60 

72367 . 1 83 7 1 M/ 1 -671 7. *2*9 31 64 

^21. 95825434 

732.9242)770 

-92.84171057 

415.70947206 

-152.31261474 

78.66362500 

24.51076870 


3.80 Diverged 


a Theoretical value of L # is one. 



TABLE II.- NON-DIMENSIONAL ANGULAR MOMENTUM AS A FUNCTION OF TIME (EXAMPLE 2) FROM INSERTION INTO TRANSLUNAR ORBIT FROM AN EARTH PARKING ORBIT 

[a - 13 U 467 n. mi. ; e = 0 . 9375 ; t Q - T = 0. 03941 hr; p = 0.2280 hr] 


Time, hr 

a L # (6 terms) 

L ' (l8 terms) 

T.' (30 terms) 

L * (k 2 terms) 

L* (54 terms) 

l' (66 terms) 

L* (78 terms) 

L # (90 terms) 

L* (102 terms) 

0.00 

1 .oanooouf) 

l.oonooooo 

1 .*0000000 

I.OOUOOOOO 

1.00000000 

1.00000000 

1,00000000 

1.00000000 

1.00000000 

0.01 

1 .00000003 

Uuooooono 

1 .00000000 

1 .00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

0.02 

l .000001112 

1.00000000 

i. 00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

0.03 

l ,000007*3 

l.uOOOOOOO 

1 .00000000 

1 .nonooooo 

1.00000000 

1.00Q0Q00Q 

UOQQQQOGQ 

WaQQGOQGQ 

1.00000000 

0.04 

t. 00003169 

1,00000000 

1.00000000 

1 .ooonoooo 

1.00000000 

1.00000000 

1.00000000 

l.oonooooo 

1.00000000 

0.05 

1.00009536 

l.oonooooo 

l.ooooonoo 

1.00000000 

1.00000000 

1,00000000 

1,00000000 

1.00000000 

1.00000000 

0 

0 

1.000233*0 

l.'JOOOOOOO 

1.000090*0 

1.00000000 

1.00000000 

1.00000000 

1,00000000 

1.00000000 

1.00000000 

0.07 

1 .00049819 

l.uOOOOOOO 

1 .*0000000 

1.00000000 

1.00000000 

1 .00000000 

1.00000000 

1 .00000000 

1.00000000 

cc 

0 

0 

1.00096694 

1 .O0O00000 

) .00000000 

1.00000000 

l.oooonooo 

1.00000000 

l.nooooooo 

1.00000000 

1.00000000 

0.09 

1.00169867 

I.oOmQOOOI 

1 .00000000 

1 .00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

0.10 

1.00283279 

1 «o0*»0O006 

1.00000000 

l.oooooono 

1.00000000 

l.nooooooo 

1.00000000 

1.00000000 

1.00000000 

0. 11 

1 .00449202 

1 .v)Q"00032 

1 .00000000 

1,00000000 

l.oooonooo 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

0.12 

1.00683246 

1.00*0014] 

l.OOOOUOOft 

l.oooonooo 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

0.13 

1 ,r>lf)03So? 

1 . o 0'<00562 

0.99999999 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

) .00000000 

l.nooooooo 

0.14 

1 .01430614 

1.U0O02013 

0.99999991 

l.nooooooo 

1 .ooonoooo 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

0.15 

1 .0198/834 

1 . 0 0006609 

0.99999933 

1,00 ft 00000 

l.nooooooo 

1.00000000 

1.00000000 

l.oonooooo 

1.00000000 

0.16 

1 .n?7010»8 

l.U0'*20097 

0.99999563 

1 .00000002 

1 .oooonooo 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

0.17 

1 ,03599024 

1.00^57129 

0.99997476 

1,00000029 

1 .0000(1001 

1,00000000 

l.nooooooo 

1.00000000 

1.00000000 

0.18 

1.04713077 

1.00152979 

0.99966804 

1. 00000294 

1.000000)9 

0.99999998 

l.nooooooo 

l.ooooonoo 

1.00000000 

0.19 

1 .06077532 

l.nQl88382 

0.99936927 

1.00002636 

1,00000332 

0.99999941 

1 .00000003 

1.00000000 

1.00000000 

0.20 

1.0772-J61 1 

1.00*39948 

0.90721876 

1.00021 091 

1.00005083 

0.99998349 

1,00000127 

1.00000036 

0.99999988 

0.21 

1.0970958) 

1.02178656 

0.98859564 

1.00152308 

1.00068187 

0.99960797 

1,00005312 

1.00002833 

0.99998387 

0.22 

1.12060894 

i. 04-55794 

0.85622234 

1,01002099 

1.00910513 

0.99197037 

1.00186145 

1.00179572 

0.99823757 

0.23 

1 .)463o3H7 

1,10443420 

0.84174949 

1,06057271 

1.08629204 

0.85623690 

1.05562255 

1.09464784 

0.84371270 

0.24 

1 . 1 80684 H 3 

1 .21 /41223 

0,45857647 

1,33887684 

1.83074814 

-1.27562261 

2.43728604 

5.21477953 

-10.44409756 

0.25 

1 .?]H29?66 

1.43*33618 

-0 . 7607560 7 

2.76759534 

8.30096441 

-31.11421387 

33,82119273 

163.75213523 

-689.15565288 


H 

-3 


Theoretical value of L is one. 



TABLE III. - NON-DIMENSIONAL ANGULAR MOMENTUM AS A JUNCTION OF TIME (EXAMPLE 3 ) FOR INITIAL CONDITIONS 
TAKEN APPROXIMATELY 46 HOURS FROM INSERTION INTO TRANSLUNAR ORBIT 

[a = 133 311 n. mi. ; e = 0 . 9762 $ t Q - T = 45-71 hr; p = 45-71 hr] 


Time, hr 

0,00 

& L * (6 terms) 
1 . OOOOOOnn 

L # (18 terms) 
1 .uO 100000 

L # (30 terms) 
1 ,00000000 

L # (4-2 terms) 
1 .non nooo n 

2.0 0 

0.99994946 

A .no-100000 

1.00000000 

l .nooooooo* 

4. 00 

0.999993*5 

l.oooonono 

1 .oonoooon 

l . 0 0 0 n n 0 0 0 

6.00 

0.99999133 

1 .00**00000 

l .nooooooo 

1. nooooooo 

8,00 

0 .9999*366 

1.00 100000 

1 .oonononn 

t .oonononn 

10.00 

0.99988968 

1 • 00 000000 

1 .nooooooo 

1 .nooooooo 

12.00 

0.99972691 

I.O0000000 

l , no non 000 

l.oonooooo 

14.00 

0.99941283 

i.oo*»onooo 

1 .nooooooo 

1 ,00t)00000 

16.00 

0,99886116 

1 .1*0900000 

l.oooonono 

l .ooonnono 

16,00 

0 .99795840 

It . 99499997 

1 .oooOoono 

l .oooonooo 

20.00 

0, 9965603* 

n. 99999985 

1 .nooooooo 

1 .00000000 

22.00 

0.9944*900 

u. 99^99924 

1 .nooooooo 

1.00000000 

24.00 

0.99152846 

0. 99999670 

1 .nooooooo 

l .nooooooo 

26.00 

0.987 424 11 

0.49498724 

0,99999999 

1 .nooooooo 

28.00 

0.981877^8 

0,99495545 

0.99949988 

l.oonooooo 

30.00 

0,97454288 

0.99485738 

0.99999915 

0.99499999 

32.00 

0.96502667 

0.99457668 

0.499994*5 

0.99999993 

34,00 

0,4528*378 

0.99-82432 

0.99996873 

0.99999914 

0 

0 

* 

vO 

r> 

0.937614*4 

o.99f*9?144 

0.99983760 

0.99990118 

38.00 

0.918663*3 

0.99235083 

0.99922893 

n.99491990 

40,00 

0 • 8954 1576 

0.48186867 

0.99662169 

0.99435044 

42.00 

0.56719430 

0.96*81007 

0.98623365 

(j. 99525338 

44.00 

0.83325970 

u ,90496386 

0,94747330 

n.96836997 

46.00 

0.792806*6 

O.H0498129 

0.81123368 

0.8063^705 

48.00 

0.74496192 

O.hl 164566 

0.35780261 

-0.09721 476 

50.00 

0,6887m3'H 

0.22 <35975 

-1.07757166 

-4.78963840 


a 


L* (54 terms) 

L* (66 terms) 

L* (78 terms) 

L (90 terms) 

L (102 terms) 

l .0000 non 0 

1. nooooooo 

1,00000000 

1 .nooooooo 

1,00000000 

1. oooonooo 

1.00000000 

l .nooooooo 

1 .nooooooo 

1.00000000 

l.oooonono 

1.00000000 

1 .00000000 

1 .oonooooo 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1 .nooooooo 

1,00000000 

1. nooooooo 

1 .00000000 

1.00000000 

1 .noonoooo 

1. nooooooo 

1 .oooonooo 

1.00000000 

1.00000000 

l.oonooooo 

1,00000000 

1.00000000 

1.00000000 

1.00000000 

1. nooooooo 

1,00000000 

1,00000000 

1.00000000 

1 .00000000 

1.00000000 

1 .nooooooo 

1.00000000 

1.00000000 

1 .nooooooo 

1.00000000 

1.00000000 

1,00000000 

1 .00000000 

1.00000000 

1 ,00000000 

1.00000000 

1 .00000000 

1 .nooooooo 

1.00000000 

1.00000000 

1 .00000000 

1 .nooooooo 

1 .00000000 

1.00000000 

1 .oooonooo 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

l. nooooooo 

1. nooooooo 

1.00000000 

1,00000000 

1.00000000 

1 .oooonooo 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

i .oooonooo 

1. nooooooo 

1.00000000 

l.oonooooo 

1.00000000 

1. oooonooo 

1 .nooooooo 

1,00000000 

1.00000000 

i.oooooooo 

0.99994998 

1 .nooooooo 

1.00000000 

1 .00000000 

1.00000000 

0.99999951 

0,49999997 

1.00000000 

1.00000000 

1.00000000 

0.99994156 

0,99999910 

0,99999990 

0.99999999 

i.oooooooo 

0,99997349 

0.99997515 

0.99999510 

0.99999903 

0.99999981 

0 . 99833935 

0,99941445 

0.99979260 

0.99992635 

0.99997381 

0.98066990 

0.90809280 

0,99263181 

0.99542841 

0.99715928 

0.79833747 

0.70830374 

0.77674350 

0.76391472 

0.74996452 

-0.90327617 

-2.32846288 

-4.84811209 

-9.30337182 

-17.18123839 

15.38392773 

-45.74771256 

133.01677686 

-384,27132685. 

1108.34246802 


m 


l Theoretical value of L is one. 



s -88 NASA- Langley, 1966 


TABLE IV.- NON-DIMENSIONAL ANGULAR MOMENTUM AS A FUNCTION OF TIME (EXAMPLE 4) FOR INITIAL CONDITIONS TAKEN AT LUNAR SPHERE OF INFLUENCE (GOING TOWARDS MOON) 


[a = -2059.46 n. mi. ; e = 1.658; t Q - T = -15-51 hr; P H = 15-51 hr] 


Time, hr 

a L # (6 terms) 

L* (l8 terms) 

L # (50 terms) 

L* (42 terms) 

L* (5^ terms) 

h* (66 terms) 

L # (78 terms) 

L' (90 terms) 

It' (102 terms) 

0.00 

1 .OOOOnono 

l.oO-'OOOOO 

1.00000000 

1.00000000 

1.00000000 

1,00000000 

1.00000000 

1.00000000 

l.nooooooo 

1.00 

1 .OO00U016 

1 .00'»00000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1,00000000 

1.00000000 

2.00 

1 .00000524 

] .00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

l.nooooooo 

3.00 

1 . O00039«0 

l.oooooono 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1 .00000000 

1.00000000 

4.00 

1 .0001*7*9 

i.oooooooo 

1.00000000 

1.00000000 

1.00000000 

1,00000000 

1.00000000 

1.00000000 

1.00000000 

5.00 

1.00051209 

1 .00 J 00 00 1 

1.00000000 

1.00000000 

1,00000000 

1,00000000 

1.00000000 

1.00000000 

1.00000000 

6.00 

1 .no 127761 

1 .U0000019 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

7.00 

1 .nn27*4*9 

1 . 0 0000260 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

8.00 

1.00539607 

1 , 00')02533 

1.00000006 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

9,00 

1.009735^4 

1 .‘J0O18Q94 

1.00000191 

1.00000002 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

1.00000000 

10.00 

1.01650765 

1. ('0114254 

1.00004114 

1.00000125 

1.00000003 

1.00000000 

1.00000000 

1 .00000000 

1.00000000 

11.00 

1.026619*2 

1 .00-383227 

1.00066446 

1.00006402 

1.00000558 

1.00000045 

1.00000003 

1.00000000 

1.00000000 

12.00 

1 .041163^1 

1. o2-»90279 

1.00848225 

1.00234740 

1.00058814 

1.00013568 

1.00002894 

1.00000560 

1.00000094 

13.00 

1,06153592 

1.10244032 

1.08918116 

1.06574330 

1,04396422 

1.02719069 

1.01550661 

1.00819230 

1.00382623 

14.00 

1.08926206 

1.36744593 

1.8001178? 

2.49128756 

3.55409557 

5.11645155 

7,30709565 

in, 23325536 

13.93487235 

15.00 

1.12621709 

2,2 1353360 

7.33652369 

30,51409329 

136.40132920 

653,94844055 

3496.4009399421382.31518555- 


16.00 

1.17454914 

4.73764381 

47.01826665 

576.36602402 

9174.39428711 

Diverged 




H 

VO 


a Theoretical value of L # is one. 
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